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Abstract — A novel channel representation for a two-hop de- 
centralized wireless relay network (DWRN) is proposed, where 
the relays operate in a completely distributive fashion. The 
modeling paradigm applies an analogous approach to the de- 
scription method for a double-directional multipath propagation 
channel, and takes into account the finite system spatial resolution 
and the extended relay listening/transmitting time. Specifically, 
the double-directional information azimuth spectrum (IAS) is 
formulated to provide a compact representation of information 
flows in a DWRN. The proposed channel representation is 
then analyzed from a geometrically-based statistical modeling 
perspective. Finally, we look into the problem of relay network 
tomography (RNT), which solves an inverse problem to infer the 
internal structure of a DWRN by using the instantaneous double- 
directional IAS recorded at multiple measuring nodes exterior 
to the relay region. 

Index Terms — Wireless relay networks, decentralized relays, 
information flows, cooperative communications, information az- 
imuth spectrum, relay network tomography 

I. Introduction 

Recently a great deal of research has been devoted to 
cooperative wireless communications Hl-Ill, which reap spa- 
tial diversity benefits from a virtual antenna array formed 
with multiple relay nodes. In a cooperative scheme, one user 
terminal (UT) partners with another UT to send its signal to the 
access point (AP) or some other final destinations. The partner 
UT serves as a relay, forwarding the message from the source 
to the destination. For most of the existing works, implemen- 
tation of the schemes assumes central coordination among 
relays. Nevertheless, the overhead to set up the cooperation 
may drastically reduce the useful throughput |4|. Furthermore, 
a centralized operation may not be feasible for a distributed 
system with a dynamic infrastructure (e.g., a wireless ad 
hoc network) or with a energy-constrained operational mode 
(e.g., a wireless sensor network). In this paper, we consider 
a decentralized wireless relay network (DWRN) without link 
combining and joint scheduling among relays |5|. 

In our previous work |6|, the performance of a DWRN 
was analyzed by realizing the noteworthy analogy between 
a virtual DWRN and a physical propagation channel. The 
information delivered from the source to the destination flows 



through a DWRN just as the signal power flows through a 
physical channel. The concept of information azimuth-delay 
spectrum (IADS) was defined, which is parallel to the power 
azimuth-delay spectrum (RA.DS) used in describing a single- 
directional physical channel fT), H). However, the method- 
ology in |6| is only applicable to a perfect receiver, which 
has an infinite angle-resolving capability and incurs zero data 
processing delay. Furthermore, the relay is supposed to listen 
and forward the data within very short intervals, comparable 
to the propagation delays of wireless signals in a DWRN. The 
current work revisits the analysis in |6| by looking into more 
general operation scenarios with finite system sensitivity and 
extended relay listening and transmitting periods. This can be 
achieved by using a double-directional description of the relay 
network, parallel to the double-directional characterization of 
the propagation channel |8|. Moreover, the DWRN is sampled 
at uniformly-spaced angles-of-departure (AODs) and angles- 
of-arrival (AOAs), leading to a discrete information azimuth 
spectrum (IAS). The problem of interest is thus to determine 
the parameters of the discrete IAS to achieve both accuracy 
with respect to the continuous benchmark model and coher- 
ence with respect to the resolution limit of the systems. Sub- 
sequently, the novel channel representation will be analyzed 
from a geometrically-based statistical modeling perspective 
El, which assumes certain geometric distributions of relay 
nodes and then derives the IAS by applying the fundamental 
information theory. 

Provided with the network description preliminaries, the 
concept of relay network tomography (RNT) is introduced, 
which is applied to identify relay locations in an unknown 
DWRN. Probing signals are transmitted from several exterior 
sources to illuminate the service region, and relayed signals 
at several exterior destinations are measured. As a result, 
the spatial distribution of active relays is obtained from the 
received instantaneous IAS by solving an inverse problem, 
which provides critical insight into the performance bottle- 
necks in a network. Apparently, the basic principle of RNT 
bears a strong resemblance to other inverse problems, in which 
key aspects of a system are not directly observable. The 



term tomography is coined to link the problem of interest 
here, in concept, to other processes that infer the internal 
characteristics of an object from external observation, as done 
in medical tomography ITOll . 

The rest of the paper is organized as follows. In Section 
II, we introduce the system model of a DWRN with realistic 
operation scenarios and formulate several principal quantities 
for description of such networks. Subsequently, the proposed 
channel representation under the framework of geometric 
network modeling is investigated in Section III, for which 
the relevant channel quantities are derived. In Section IV, we 
provide a general formulation of the RNT problem. Section V 
demonstrates the properties of the model parameters and the 
efficacy of the RNT methodology through several numerical 
examples. Finally, some concluding remarks are drawn in 
Section VI. 

II. Double-Directional Description of DWRN 

A. System Model 

We will consider a DWRN with L + 2 nodes: A source S, a 
destination D, and L non-collocated relays Ri, R2, • • • , Rl 
as illustrated in Fig. 1. The relays are randomly scattered in 
a given region TZ, which are due to the irregular nature of 
wireless sensor or ad hoc networks. In a sensor network, 
a large number of relays with sensing capabilities may be 
spread over the site under investigation in a distributive and 
randomized manner. On the other hand, each UT in an ad 
hoc network can serve as an assisting relay and the non- 
stationary network architecture may cause uncertainty to the 
relay positions 1, 1 1 J . Subsequently, the following assumptions 
on the transmission strategy are imposed. 

1) Relay Distribution: The locations of R; {I = 
1, 2, • • • ,L) are described by a statistical density function 
/r (R-i) (R-i G T^), which represents the ensemble of a 
large number of randomly distributed relays. 

2) Relay Selection and Cooperation: The network geometry 
is unknown to either S or D. Thus, S broadcasts the same 
message to all R; and coherent relaying is not employed 
(i.e., without link combining and joint scheduling among 
relays). 

3) Decode-and-forward Processing: Each R; receives infor- 
mation from S in the first hop (backward channel) and 
then forwards the decoded signal to D in the second hop 
(forward channel). Two orthogonal frequency tones are 
used at these two transmission phases, respectively. 

4) Medium Access: Due to the decentralized and unregu- 
lated nature of relay positions, a simplified time-division 
protocol is applied, which requires no knowledge of the 
network geometry at either S or D. At the initial time 
t = to, S transmits its signal to all R; = 1, 2, • • • , L). 
For the ^th relay, it listens to the transmission from S 
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Fig. I. Double-directional DWRN channel 



in the backward and forward phases, respectively. The 
same time-division process repeats in the subsequent 



time slots 



At 



(B) 
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to + kx (AtP' + Atf^^ , to + (fc + 1) X 



(fc = 1, 2, • • • ). For simplicity, each 
relay spends equal times in its listening and forwarding 
phases. Apparently, R( introduces an AOD fi; and an 
AOA during each transmission cycle as shown in Fig. 
1 . It is worth emphasizing that the time-of-arrival (TOA) 
Ti — ^'^'^^ ^ which is one of the key quantities 
used in the multipath channels [8J, has lost its physical 
meaning in the current context. Here ds,Ri and c/r, ,d 
are the distances of the links S — > R; and R; — > D, 
respectively, and c is the speed of electromagnetic waves. 
5) Frequency-Tone Assignment: In general, mutual interfer- 
ence among the forward channels may incur due to the 
simultaneously transmitting relays. One way to eliminate 
the interference is to assign L frequency tones drawn 
from the pool of available bandwidth [w„iin, Wmax] to L 
relay-to-destination pairs, all of which are also orthogonal 
to the carrier allocated to the backward channels. 

B. Double-Directional Representation of DWRN 

For a two-hop DWRN, multiple information pipelines are 
laid via the L Hnks S Ri -> D, • • • , S Rl D. 
Furthermore, the intersection of each AOD and each AOA 
uniquely determines the location of a relay. Therefore, as an 
ideal reference model, the information flows in the network 
are continuously distributed in the AOD-AOA domain as 



(1) 



(=1 



parameters AtJ ' and AtJ ' denote the time spent 



where /D(r2,^) is defined as the double-directional IAS and 
/; is the outage capacity given an outage probability of Pout 
for the link S — > R/ — > D. The terminology double-directional 
is adopted according to |8|. 

Eq. (1) assumes perfect receivers that are able to resolve 
the departing and arriving information flows with infinite 



sensitivity in the angular domain. Therefore, the formula- 
tion represents an ideal benchmark model. In practice, both 
the source and destination have limited system sensitivity. 
Consider two data streams via two relays R;j and R/^ 
(l<^i7^^2<i)- Due to the finite length of the transmit (receive) 
antenna aperture, the AODs (AOAs) of these two information 
flows cannot be successfully distinguished if the separation of 
their departing (impinging) angles, iJl/j-Sl/J (l^f/^ is 
within the resolution limit of the source (destination) antenna, 
AQ (A^f). Therefore, die DWRN channel should be sampled 
at uniformly-spaced AODs and AOAs: 
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where [-J and [•] are the floor and ceiling functions, respec- 
tively. 

Subsequently, the discrete double-directional IAS can be 
expressed as 

{n,, VI/,) = /D(r!, ^)fl;;i\n, ^)dnd^, 
(17,*) : n 



where = 



In (4), The superscript 



(d) 



An ^ ^ An 
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(4) 

denotes a discrete model. 
The local joint AOD-AOA pdf f^;l\n,^) satisfies 

, f^'{,\n,'i')dnd"i' — 1, where the integration domain 
Si J- is illustrated in Fig. 1. 

III. Geometrically-Based Statistical Model for a 
Random DWRN 

We will discuss the proposed analytical framework in 
Section II from a geometrically-based statistical modeling 
perspective [6J, where a large number of relays are ran- 
domly located in the two-dimensional space according to 
a specified relay density function. This approach is useful 
when the spatial structure observed in a large DWRN is 
far from being regular ifTTl . The medium access arrange- 
ment follows the protocol in Section II-A. Moreover, all 
channels experience independent frequency-flat block fading 
with the complex channel gain between the link M — > N, 
/iM,N (Me {S,Ri,--- ,RL},Ne {D,Ri,--- ,Rl}) fol- 
lowing the Nakagami distribution with shape parameter m 
JS). The associated instantaneous power |/im,nP is gamma- 
distributed with the same shape parameter, i.e., |/im,nP ^ 
, where Q{k,9) represents the gamma distribu- 



tion with scale 9 and shape n. The mean value Am,n = n 
with ly being the path loss exponent. The ambient noise at the 
relays and the destination, ^ A/'(0,iVo), where Af{^,a'^) 
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Fig. 2. Geometrically-based statistical model for a DWRN 



represents a normal distribution with mean /i and variance 
a^. The communication bandwidth is B for any channel. 
Subsequently, the average received signal-to-noise ratio (SNR) 
for each link is SNR = -g^, where P is the average transmit 
power for all transmitting terminals. Finally, the channel state 
information is assumed to be known at the receivers but 
unavailable to the transmitters Ol. 

Let the x-y coordinate system be defined such that the 
destination is at the origin and the source lies on the x axis, 
as shown in Fig. 2. It is assumed that there are many relay 
nodes, the locations of which are described by the statistical 
relay density function fx.y{x, y). By applying the law of sines 
to the triangle SRD in Fig. 2, the relay-to-destination distance 
c^R.D can be expressed in terms of the AOD n and the AOA 
^ as 

ds.D sinfi 



'^^■°-sin(l] + vi/) 
The source-to-relay distance c?s,r is similarly derived as 

ds.D sin ^E* 



sin(r2 + *) 



(5) 
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Proposition 1 (Derivation o/ /j2 ,f (17, 4')).- 
The joint AOD-AOA pdf of information flows within a 
DWRN is given by 



/o,*(f^,*) 



|sin2fi + sin2^'-sin2(fi + *)| 

[1 - 008 2(17 + *)]^ 

ds D sin il COS vl* ds.D sin il sin 



sm{n + *) ' sm{n + 
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where ds,r> is the distance between the source and the 
destination. 

Proof: The proof is omitted here for simplicity. 

Proposition 2 (Derivation o//d(17, vf)).' 

The double-directional IAS Id (17, vp) is the argument of the 



function 
V{I) 



7 I m 
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r(m) 



(8) 
is the 



when V{I) = Pout. In (8), 7(a,/3) 
lower incomplete gamma function and r( ) is the gamma 
function. 

Proof: The proof is omitted here for simphcity. 
Finally, the local joint AOD-AOA pdf in (4), f^l^{^, 
is related to /n,*(0, ^) through the following relationship: 



(9) 



Substituting (7)-(9) into (4) gives the discrete double- 
directional IAS. 

IV. RNTforaDWRN 

A. Formulation of RNT 

In the general RNT problem depicted in Fig. 3, L relays 
(Ri,-- - ,Rl) are randomly scattered in a finite area TZ 
and Q measuring nodes (Ni,-- - ,Nq) are placed exterior 
to TZ. Nevertheless, the following analysis is also appli- 
cable to any other measurement network orientation with 
respect to TZ. Each measuring node can serve as the probing 
source as well as the information destination, which has a 
transmit/receive azimuthal resolution of AQ. As shown in 
Fig. 3, the first node Ni transmits over its full spectrum 
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denotes Ni. All other nodes 



N2 



Nq record an information flow. Each receive 
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superscript denotes Ng. At each jAQ, Ng 
whole frequency band within which the relays are operating. 
When Ng locks on the Ith. frequency component wi, which 
yields a detectable information flow associated with the relay 
R;, it then monitors the instantaneous channel capacity over an 
observation period and records the time- variant capacity I (t). 
From the set of 7 {t), the channel outage capacity can 
be estimated empiricaUy, where the superscript ^^'^^ denotes 
the information pathway from Ni to N^. The similar process 
continues until all the active relays are identified. Then the 
second node N2 sends the probing message and all the other 
nodes Ni, N3, • • • , Ng receive the data, and so on. 

A procedure to calculate the double-directional IAS for a 
DWRN with given sources, relays and destinations is called 
the direct problem. A procedure to obtain relay distribution 
knowing sources/destinations and measuring the information 
flows at multiple receivers is called the inverse problem. The 
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Fig. 3. Formulation of RNT 



direct problem has already been discussed in Section III. The 
corresponding inverse problem can be formulated as follows. 

RNT Problem: Given the set of measuring nodes 
Ni,--- ,Nq, the estimated (Q — l)^xL single-directional 
AOA matrix 



* = [ *i 



*2 

(1,2) 



(1.2) 



(1,2) 



(Q,Q-i) 



(Q,Q-i) 



(10) 



and the estimated {Q — l)^xL outage capacity matrix 



[ii 



f(l,2) 



f(l,2) 



jiQ.Q^l) f(Q,Q-l) 
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f(l,2) 
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(11) 
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and /, 



(91,92) 



are the empirical AOA and channel 



outage capacity associated with the information pipeUne 



N, 



91 



R; 



Nq2, respectively. Note that due to the 



channel reciprocity in a two-hop DWRN, it is expected that 

r(92,9l) 



= n^'^'"^' and il""''''^ = II"'-"". Therefore, (11) is 



equivalent to a double-directional AOD-AOA matrix, where 
half of the AOAs correspond to the AODs of the reversed 
information pathways. Identify the locations of Ri, • • • ,Ri 
that satisfy the foUowing two constraint functions 



ei 



and 



<ei, 

! (i; - II) < €2, 



1 = 1,2,- 



1 = 1,2, 



(12) 
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and 1; are the AOA and channel outage capacity vectors 
solved in the direct problem for all the information paths via 
R; assuming a specific location of R;. 



B. Algorithms to Solve The RNT Problem 



Hk^, which is given by 



In the current work, we will 
norm of the two error vectors, i.e.. 



consider 



the 



E,.E,.^,. (^(..^..)_^(..^.^)) and 



62 \ll - 1/ j = 

The objectives are to 



minimize the values of 62 and meanwhile ensure that ei = 
for all the detectable information flows. The optimization 
process can be reaUzed through the following steps: 

1) Discretize the solution domain of the relay locations 
TZ into W sufficiently small square cells, i.e., R = TZ^'^^ = 
{xi,--- ,yiw}, where y.^ = {x^,y^) {w = 1,2,--- ,W) 
represents the Cartesian coordinates of the center of the wth 
cell; 

2) For each identified information pipeline N^^ ^ R; ^■ 



, determine the subset of '. 



fiqim) ,1,(91,92) 



(91,92) 



^(91,92) 



C R, such that Vx e 



(x); 



3) Determine the joint of all the sets obtained in Step 2 as 

ll = Rp'2) pi ]j(l,3) pi . . . pi m(Q'Q-I) ; 

4) Estimate the location of R; as 



x(Ri 



arg mm 
xeR( 




(91,92) 7-(9l,92) 



'(x) 



(14) 

It is not always possible to obtain an accurate estimate of 



the outage capacity /, 



(91,92) 



due to the limited observation 



data at each destination. In such a case, the inverse problem 
could be approached from a hypothesis-testing perspective. An 
alternative optimization procedure to Step 4 can be formulated 
as follows. 

4*) Let there be K relay-position hypotheses denoted 
Hi : x(R;) = xi,'H2 : x(R;) = X2, • ■ ■ Mk ■ x(R;) = x^, 
where K = and R; = {xi,X2,--- ,^k}- Furthermore, 
let eki,k2 be the desired probability of incorrectly selecting 
Hk2 given that Hki is true. The sequential probability ratio 
test procedures will be applied for statistical decision making, 
where the decision at each stage is based on the likelihood 
ratio. In particular, the multi-hypothesis sequential test is 
considered and summarized below |fT2| . The likelihood ratio 
including prior probabilities for a pair of hypotheses ki and 
^2 after the oth data observation is 
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/«.,(Ji)/?^.,p2)---/«.,(Jo)Pr(HfcJ 
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where (a) follows the fact that the instantaneous outage 
capacities are sequentially observed for all the information 
paths via R; during the measurement process. Therefore, they 
can be assumed to be independently distributed. The equality 
in (6) is obtained by differentiating the cdf in (8) with respect 
to the variable /. In (16), 



{it 



TO 4 



1 



and 



„ / 7-(9l,92) -- 



SNRx [ds,R(xfcJ]" 



TO f 4^',° - Ij 
SNRx [dR^D (xfej]" 



If the number of observations is not fixed, a decision on 
the relay-position state (i.e., T-L-^. is selected) is made when the 
condition 

1 fMl VfcT^fca (17) 



-k,k2 



is met. The threshold in (17) to which the likelihood ratio 
is compared is designed to ensure that the average rate of 
making a wrong decision in favor of "H^j when Hki is true is 
no more than eki,k2- Nevertheless, the inequality in (17) may 
not always be satisfied due to the insufficient observation data. 
In such a case, a hypotheses is selected if its associated 
likelihood has the maximum value after a fixed number of O 
observations, i.e.. 



where 3o 



/^(1.2) ^{1, 



3) 



r(Q,Q-l) 



is the instanta- 



neous outage capacity vector at the oth observation for all the 
information routes via R;. f-u^ (3o) is the pdf of 3o under 



fc = argmax (3i)/« J32)- 



•/«,po)Pr(Hfe) (18) 



Eqs. (17) and (18) provide alternative objective functions to 
the one given in (14). 
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Fig. 4. Double-directional IAS: (a) continuous reference model and (b) 
discrete model for systems with an angular resolution of 10° 
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Fig. 5. Estimation of relay locations using the proposed RNT algorithm 



V. Numerical Examples 

We consider the high SNR regime (SNR = 30 dB) for a 
fixed outage probability Pout = 0.01. The shape parameter 
of the Nakagami distribution is m = 1 (Rayleigh fading). 
The distance between the source and the destination rfs.D 
is assumed to be 100a/3 m (see also Fig. 2) and the path 
loss exponent is v = —3. The relay region is assumed to 
be a circle centered at O(50\/3 m, 50 m) with a radius 
r = 40 m. Finally, the spatial resolutions of both the source 
and destination antennas are 10°. 

Fig. 4 illustrates the double-directional IAS for both the 
continuous reference model in (1) and its discrete counterpart 
in (4). Apparently, both Ir>{il,'i>) and ll^\ft,'^) decrease 
as the AOD O and the AOA increase. This certainly 
holds an intuitive appeal because larger values of Q and 4* 
lead to more severe path loss, thereby reducing the outage 
capacity of a relay path. As a practical concern, Fig. 4 
provides an intuitive and quantitative answer to the following 
two important questions: how will the information flows be 
scattered in the AOD and AOA domains and what will be the 
beam-scanning ranges of the source and destination antennas 



in order to achieve sufficiently high outage capacities? 

Next, we look into the RNT problem discussed in Section 
IV. A three-node measurement network is deployed as illus- 
trated in Fig. 5. A number of relays are uniformly distributed 
within a circular region. The procedures presented in Section 
IV-B are applied to estimate the relay locations, where the 
objective function of (18) is used. The number of observations 
of the instantaneous channel capacity for each information 
pipeline is O = 10. Furthermore, no a priori knowledge of the 
probability of each location under hypothesis test is available. 
As can be observed from Fig. 5, most of the estimated relay 
locations agree well with the actual ones, which demonstrates 
the efficacy of the proposed RNT algorithm by using a 
relatively small number of measuring nodes and observations. 

VI. Conclusions 

We have presented a double-directional description method 
for analyzing two-hop DWRN channels, which is based on 
an analogous approach to the double-directional physical 
wireless channels. The proposed model has been investigated 
using a geometrically-based statistical modeling approach with 
uniformly distributed relay nodes. We have also investigated 
the RNT problem by applying the proposed network rep- 
resentation. Finally, numerical examples have been used to 
demonstrate the practical significance of the proposed analyt- 
ical framework. 
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